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O ■ Abstract 
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Using the two-loop approximation of chiral perturbation theory, we calculate the mo- 
mentum and density dependent nuclear spin-orbit strength Ui s (p, kf). This quantity is de- 
rived from the spin-dependent part of the interaction energy T, sp i n = § <?• (<f xp) Ui s (p, kf) 
of a nucleon scattering off weakly inhomogeneous isospin symmetric nuclear matter. We 
CsJ | find that iterated l7r-exchange generates at saturation density, kfo = 272.7 MeV, a spin- 

orbit strength at p = of t/; s (0, kfo) ~ 35 MeVfm 2 in perfect agreement with the empirical 
value used in the shell model. This novel spin-orbit strength is neither of relativistic nor 
^sD ! of short range origin. The potential Vi s underlying the empirical spin-orbit strength 

Uis = Vis r f s becomes a rather weak one, Vi s ~ 17 MeV, after the identification ri s = m^ 1 
as suggested by the present calculation. We observe however a strong p-dependence of 
Ui s (p,kfo) leading even to a sign change above p = 200 MeV. This and other features of 
the emerging spin-orbit Hamiltonian which go beyond the usual shell model parametriza- 
tion leave questions about the ultimate relevance of the spin-orbit interaction generated 
by 27r-exchange for a finite nucleus. We also calculate the complex-valued isovector single- 
particle potential Ui(p, kf)+i Wi(p, kf) in isospin asymmetric nuclear matter proportional 
to Ts(N — Z)/(N + Z). For the real part we find reasonable agreement with empirical 
■ values and the imaginary part vanishes at the Fermi-surface p = kf. 
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1 Introduction and summary 

The introduction of the spin-orbit term into the nuclear single-particle Hamiltonian by Haxel, 
Jensen, Suess and Goeppert-Mayer [|TJ in 1949 has been most decisive for the success of the nu- 
clear shell model. Only with a very strong and attractive spin-orbit potential one is, for example, 
able to explain the observed sequence of so-called magic numbers {2, 8, 20, 28, 50, 82, 126, . . .}. 
The dynamical origin of the strong nuclear spin-orbit force has not been fully resolved even up 
to date. The analogy with the spin-orbit interaction in atomic physics gave the hint that it could 
be a relativistic effect. This idea has lead to the construction of the relativistic (scalar- vector) 
mean- field model 0. In this model the spin-independent nuclear potential (of approximate 
depth —50 MeV) results from an almost complete cancelation of a very strong attraction gen- 
erated by scalar (u-meson) exchange and a nearly equally strong repulsion generated by vector 
(cj-meson) exchange. The corresponding spin-orbit term (obtained by a non-relativistic reduc- 
tion of the nucleon's Dirac-Hamiltonian) comes out proportional to the coherent sum of the 
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very large scalar and vector mean-fields. In this sense, the relativistic mean-field model gives 
a simple and natural explanation of the basic features of the nuclear shell model potential. 
Refinements of relativistic mean-field models which include additional non-linear couplings of 
the scalar and vector fields or explicitly density-dependent couplings are nowadays widely and 
successfully used for nuclear structure calculations 0, [|, || . 

The nuclear spin-orbit potential arises generally as a many-body effect from the underlying 
spin-orbit term in the (free) nucleon-nucleon scattering amplitude. The calculation of the 
tree level diagrams with one scalar-meson or one vector-meson exchange between nucleons 
gives indeed a spin-orbit term in the NN T-matrix proportional to 1/M 2 , with M denoting 
the nucleon mass. The nuclear spin-orbit potential corresponding to scalar and vector meson 
exchange is therefore obviously a truly relativistic effect. However, the quadratic reciprocal 
scaling of the spin-orbit NN-amplitude with the nucleon mass M is not universal, and it changes 
if one considers the exchange of two mesons between nucleons, i.e. loop diagrams. For example, 
irreducible two-pion exchange gives rise to a spin-orbit term in the NN T-matrix proportional 
to 1/M (see eqs. (22,23) in ref. 0) and iterated one-pion exchange produces a spin-orbit term in 
the NN T-matrix which even scales linearly with the nucleon mass M (see eq.(33) in ref.||). It 
is one of the chief purposes of this paper to investigate in detail the contributions from iterated 
one-pion exchange to the nuclear spin-orbit interaction. As already mentioned the latter arises 
from the spin-orbit term in the NN T-matrix as a many-body effect, e.g. in connection with a 
nuclear matter calculation. 

In a recent work 0, we have used chiral perturbation theory for a systematic treatment 
of the nuclear matter many-body problem. In this calculation the contributions to the energy 
per particle, E(kf), originate exclusively from one- and two-pion exchange between nucleons 
and they are ordered in powers of the Fermi momentum kf (modulo functions of kf/m^). It 
has been demonstrated in ref.|7| that the empirical saturation point and the nuclear matter 
compressibility K ~ 255 MeV can be well reproduced at order 0(k 5 j-) in the chiral expansion 
with just one single momentum cut-off scale of A ~ 0.65 GeV which parametrizes all necessary 
short range dynamics. Most surprisingly, the prediction for the asymmetry energy, A(kfo) = 
33.8 MeV, is in very good agreement with its empirical value. Furthermore, as a nontrivial 
fact pure neutron matter is predicted to be unbound and the corresponding equation of state 
agrees roughly with that of sophisticated many-body calculations for low neutron densities 
Pn < 0.25 fm -3 . In a subsequent work ||, the momentum and density dependent (real) single- 
particle potential U(p,kf) (i.e. the spin-independent average nuclear mean-field) has been 
calculated in the same framework. It was found that chiral In- and 27r-exchange give rise to a 
potential depth for a nucleon at the bottom of the Fermi sea of U(0, kfo) = —53.2 MeV. This 
value is in very good agreement with the depth of the empirical optical model potential and the 
nuclear shell model potential. Nuclear matter at finite temperatures has been investigated in the 
same framework in ref.0. There it was shown that chiral In- and 27r-exchange reproduce the 
first-order liquid-gas phase transition of isospin symmetric nuclear matter with a realistic value 
T c ~ 19 MeV of the critical temperature. Our approach to the nuclear matter problem is in 
many respects different from most other commonly used ones, where one starts from a so-called 
realistic NN-potential. For example in the relativistic nuclear matter calculation of ref . |10] the 
S-, P— and D- waves deliver more than 95% of the potential energy per particle. The finding 
that perturbative chiral pion-nucleon dynamics leads already to good nuclear matter and single 
particle properties hints at the fact that the detailed NN-interaction is of no more relevance. 
Fine-tuning of the single cut-off scale A to one nuclear matter observable (the binding energy 
per particle —E(kf )) is however still necessary in our present approach 

It is the purpose of this work to calculate, using the same framework as in ref.||, the mo- 
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mentum and density dependent nuclear spin-orbit strength Ui s (p, kf). This quantity is derived 
from the spin-dependent part of the interaction energy S sp j n = | a- (qxp) Ui s (p, kf) of a nucleon 
scattering off weakly inhomogeneous isospin symmetric nuclear matter. We will present here 
analytical expressions for the contributions from l7r-exchange and iterated l7r-exchange to the 
spin-orbit strength Ui s (p,kf). Furthermore, we calculate in isospin asymmetric (homogeneous) 
nuclear matter the (complex-valued) isovector single-particle potential Ui(p,kf) + iWi(p,kf) 
accompanied by the isospin double-asymmetry r 3 (iV — Z)/(N + Z). Our results can be sum- 
marized as follows: 

i) At nuclear matter saturation density, kfo = 272.7 MeV, l7r-exchange and iterated lir- 
exchange generate a spin-orbit strength at p = of Ui s (0, fc/o) = (0.4 + 34.7) MeVfm 2 . 
This result, which is dominated by the contribution of four Hartree-type diagrams, is in 
perfect agreement with the empirical value of the spin-orbit strength U[ s ~ 35 MeVfm 2 
|TT| , P^ ] used in shell-model calculation of nuclei. The novel spin-orbit strength found here 
is neither of relativistic nor of short range origin. It is in fact linearly proportional to the 
nucleon mass M and its inherent range is the pion Compton wavelength m" 1 = 1.46 fm. 
The latter feature tempts to an unconventional interpretation of the strong nuclear 
spin-orbit interaction. The potential Vi s underlying the empirical spin-orbit strength 
Uis = Vi s rf s ~ 35 MeVfm 2 becomes a rather weak one, namely Vi s ~ 17 MeV, after the 
identification of the effective range r/ s with the pion Compton wavelength, ri s = m" 1 , as 
suggested by the present calculation. 

ii) We observe however a strong p-dependence of Ui s (p, kf ) which leads even to a sign change 
above p = 200 MeV. The calculated spin-orbit strength U[ s (0,kf ) depends also strongly 
on the value of the pion mass and it shows a pronounced maximum around m n ~ 100 MeV. 
A further property of the spin-orbit Hamiltonian emerging from our diagrammatic calcu- 
lation is that it has (in coordinate space) terms proportional to V/(r) as well as terms 
proportional to /(r) V/(r) (with f(r) = p(r)/p(0) the normalized radial density profile) 
which get weighted differently at the surface of a finite nucleus. All such features of 
our calculation which go beyond the usual shell model parametrization of the spin-orbit 
Hamiltonian leave questions about the ultimate relevance of the spin-orbit interaction gen- 
erated by 27r-exchange for a finite nucleus. Implementing the present results for Ui s (p, kf) 
into nuclear structure calculations will clarify the role of the nuclear spin-orbit interaction 
generated by 27r-exchange. 

iii) The real part of the isovector single-particle potential Ui(0,kf) generated by chiral ln- 
and 27r-exchange has a density dependence very similar to that of the asymmetry energy 
A{kf) J7J. At saturation density, p — 0.178 fm~ 3 , we find a repulsive isovector potential 
of Ui(0,kf ) = 47 MeV. This prediction is comparable to the value U\ ~ 33 MeV used 
in shell model calculations iTll or the value U\ ~ 40 MeV deduced from nucleon-nucleus 



scattering in the framework of the optical model fL3[ . The momentum dependence of 
Ui(p, kfo) is non-monotonic in the interval < p < kfo- One observes a broad maximum 
at p = 230 MeV where the (real) isovector single-particle potential has increased by about 
30% to the value 63 MeV. The imaginary part Wi(p,kf) vanishes (quadratically) at the 
Fermi-surface (p = kf) in accordance with Luttinger's theorem [|T 
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2 General considerations about the spin-orbit term 



Let us begin with recalling the spin-orbit Hamiltonian of the nuclear shell model [11] which is 
generally written in the form: 
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Here, a is the conventional Pauli spin-vector and I = f x p is the orbital angular momentum 
of a nucleon. p(r) denotes the radial density distribution of a nucleus (typically parametrized 
by a Saxon- Woods function). The empirical value of the nuclear spin-orbit strength is U\ a ~ 
35 MeVfm 2 [|TTJ, [T2|. This large and posit ive value of C7j s leads to a strongly attractive spin-orbit 
potential acting mainly at the surface of a nucleus. 

We wish to calculate the nuclear spin-orbit strength Ui s (or an appropriate generalization 
of it) in the systematic framework of chiral perturbation theory |7], ||. The first observation 
one makes is that the spin-orbit interaction vanishes identically in infinite homogeneous nuclear 
matter since there is no preferred center in this system in order to define an orbital angular 
momentum. Therefore one has to generalize the calculation of the single-particle potential in 
ref.|| to (at least) weakly inhomogeneous nuclear matter. The relevant quantity in order to 
extract the nuclear spin-orbit strength is the spin-dependent part of the interaction energy of 
a nucleon scattering off weakly inhomogeneous isospin symmetric nuclear matter from initial 
momentum p — q/2 to final momentum p + q/2, which reads: 



Z spin = - a ■ (<f x p) Ui s (p, kf) 



(2) 



The (small) momentum transfer q is provided by the Fourier-components of the inhomogeneous 
nuclear matter distribution. The density form factor $(<f) = / d 3 r e~ tq ' r f(r) plays the role of a 
probability distribution of these Fourier-components. The form factor <&(q) should be viewed as 
narrowly peaked function around q = with its Fourier-transform equal to the (slowly varying) 
density profile, f(r) = (2tc)~ 3 J d 3 qe tq r $(<f). Using this relationship the spin-orbit Hamiltonian 
Tiis in eq.(l) becomes equal to the Fourier-transform of the product of the density form factor 
and the spin-dependent interaction energy: Hi s = (2ti)~ 3 f d 3 q e iq ' r $(q) ^ sp in- Consistent with 
the assumption of a weakly inhomogeneous nuclear matter distribution we keep in T, spin only 
linear terms in q corresponding to small spatial density gradients. For practical purposes, this 
means that after isolating the proportionality factor q in an explicit calculation the momentum 
and density dependent spin-orbit strength Ut s (p,kf) can be finally computed in the limit of 
homogeneous isospin symmetric nuclear matter (characterized by its Fermi momentum kf). 
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Fig.l: Left: The double line symbolizes the medium insertion for a weakly inhomogeneous many- 
fermion system T(pj,qj) defined by eqs.(3,4)- Middle: The a- and to- exchange Hartree graph. 
Right: The lrc-exchange Fock graph. 
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In ref. f| the calculation of the single-particle potential in homogeneous isospin symmetric 
nuclear matter has been organized in the number of so-called medium insertions. The latter is 
a technical notation for the difference between the vacuum and in-medium nucleon propagator 
(see eq.(3) in ref.|7]]). In the case of homogeneous nuclear matter a medium insertion in a self- 
energy diagram converts effectively a four- dimensional loop integration into an integral over a 
Fermi-sphere of radius kf. The medium insertion for a (non-relativistic) many-fermion system 



is generally constructed from the sum over the occupied energy eigenstates as [15 



r(Pi,£) = / d 3 ri J d 3 r 2 £ ^(f^KfOe^'^^e-J*-^^. (3) 

E£occ 

The double line in the left picture of Fig. 1 symbolized this medium insertion together with the 
assignment of the out- and in-going nucleon momenta pj ± ffj/2. The momentum transfer qj is 
provided by the Fourier- components of the inhomogeneous matter distribution. Semiclassical 
expansions [TJJ give for a weakly inhomogeneous and spin-saturated many-fermion system: 

T(p 3 , q 3 ) = 6(kf - \pj\) m) {l + 0%)} , (4) 

with ~ / d 3 pjT(pj, qj) the density form factor introduced after eq.(2). The subleading 

0(qj) term in eq.(4) will in fact never come into play in our diagrammatic calculation of the 
spin-dependent interaction energy £ sp j n to linear order in q. 

As a first check on this formalism we evaluate the a- and cu-exchange Hartree diagram in 
Fig. 1. We perform the non-relativistic 1/M-expansion of the scalar/vector interaction vertex 
sandwiched between Dirac-spinors for the out- and in-going nucleon (of momentum p±q/2) 
until we obtain the spin-orbit like term iff ■ (q x p)/4M 2 . After that we can take the limit of 
homogeneous nuclear matter and perform the remaining integral over a Fermi-sphere of radius 
kf which brings one factor of density p. Putting all pieces together we reproduce the familiar 
result: 

P =% (5) 

2M l \m l a J ott z 

of the relativistic mean- field model J|, ||| . The contribution of the analogous a- and w-exchange 
Fock diagrams to the nuclear spin-orbit strength can also be easily calculated with the help 
of the formalism outlined above. We obtain from these Fock diagrams 1/4 of the cx-exchange 
contribution and 3/4 of the u;-exchange contribution written in eq.(5). In this work our main 
interest is fo cussed on the nuclear spin-orbit interaction generated by chiral one- and two-pion 
exchange. 



3 Diagrammatic calculation of the spin-orbit strength 

In this section we present analytical results for the nuclear spin-orbit strength Ui s (p, kf) as given 
by chiral one- and two-pion exchange. We start with the l7r-exchange Fock graph (last diagram 
in Fig. 1). In the static approximation the product of ^TV-interaction vertices ff ■ (p — pi) ff ■ 
(p — pi) = (p — Pi) 2 is spin-independent. A non- vanishing nuclear spin-orbit strength comes 
therefore only as a relativistic l/M 2 -correction. Isolating the iff x q factor from the product 
of fully relativistic pseudovector 7rA-interaction vertices, performing the 1/M-expansion, and 
integrating finally over a Fermi-sphere of radius kf,we get from the l7r-exchange Fock diagram: 
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L(x, u, 



(6) 



Here, we have introduced the dimensionless variables x = p/m n and u = kf/m n and the 
auxiliary function 

L(x,u ) = —\n—^ 7 

v ' Ax 1 + (u - x) 2 

For the readers convenience we will present the spin-orbit strength at p = of each individual 
diagram in a separate formula, since in most cases the limit x — > is quite non-trivial. For the 
l7r-exchange Fock diagram we find the simple expression 
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(8) 



which gives numerically at saturation density kfo = 272.7 MeV only about 1.1% of the empirical 
value Ui s ~ 35MeVfm 2 . 



(b) 



(d) 



Fig. 2: The iterated In -exchange Hartree graphs. Their isospin factor in symmetric nuclear 
matter is 6. 

Next, we come to the evaluation of the four Hartree diagrams of iterated one-pion exchange 
shown in Fig. 2. We start with the left graph with one medium insertion, labeled (a). The 
relevant id x q prefactor can be isolated already in the first step of the calculation from the 
product of 7T ^-interaction vertices d ■ (I — q/2) a ■ (I + q/2) at the open nucleon line. For all 
remaining parts of the diagram one can then take the limit of homogeneous nuclear matter (i.e. 
q = 0). Using the analytical results given in section 4.3 of ref.JJ for the inner <i 3 Z-loop integral 
we can even perform the integral over a Fermi-sphere of radius kf. Altogether, we find the 
following closed form expression for the spin-orbit strength generated by the Hartree diagram 
(a): 

jl6(w 3 + x 3 ) arctan(u + x) + 16(x 3 — u 3 ) arctan(u — x) 

xL{x, u) > . (9) 



Ut\p,k f ) 



2g\Mm\ 
+ux(7 - 9m 2 - 9x 2 



9(u 2 - x 2 ) 2 - 30m 2 - 30x 2 - 7 



At zero nucleon momentum (p = 0) this simplifies to 

g\Mm 2 n 



32tt 3 /4 



4 arctan u — 3m 
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1 + u 2 



(10) 
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Note also that l/j-g (j>,kf) in eq.(9) originates from the real part of the iterated l7r-exchange 
spin-orbit NN-amplitude (eq.(33) in ref.0) evaluated in forward direction and integrated over 
a Fermi- sphere of radius kf. 

We continue with the calculation of the Hartree diagrams with two medium insertions, 
labeled (b), (c) and (d) in Fig. 2. In these diagrams, a (small) momentum transfer gl^ 
with q = cfi + q 2 occurs at each medium insertion. The important prefactors ia x cfj can 
again be isolated in the first step of the calculation from the product of 7riV-interaction ver- 
tices at the open nucleon line. After that the vectors qj can be set to zero in all remain- 
ing components of these diagrams. When Fourier-transformed with the density form fac- 
tors to coordinate space each such momentum transfer qj , (j = 1, 2) leads to the expression: 
(27r)~ 6 / d 3 q J d 3 qj e iq ' r icfj $(<£/) $(<f — qj) = f(r) V/(r), with f(r) the density profile of weakly 
inhomogeneous nuclear matter. Consistent with the assumption of a weakly inhomogeneous 
nuclear matter (i.e. keeping only linear terms in small spatial gradients) we can make the 
replacement: f(r) V/(r) — ■> V/(r). We will come back to this point in the next section when 
discussing the results for the spin-orbit strength Ui s (p, kf) as well as their relevance for a finite 
nucleus. The essential conclusion from the previous considerations is that for the calculation 
of S sp j„ in weakly inhomogeneous nuclear matter each momentum transfer qj can be identified 
with q. With the help of this rule and certain techniques to reduce six-dimensional principal 
value integrals over the product of two Fermi-spheres of radius kf, we end up with the following 
result for the Hartree diagram (b): 

6g%Mml f 1 , y \ n , /2 2 2 ^ u + xyM3s + 2s 3 



U^( P ,kj 



with the auxiliary function s 



dy — 

1 x 



2uxy + (u — x y ) In 



u 



xy 



1 + s 2 



3 arctan s 



(11) 



xy + \Ju 2 — x 2 + x 2 y 2 . Throughout this work the momentum p 



is restricted to the interval < p < kf. At zero nucleon momentum (p = 0) eq.(ll) simplifies 
to 



U%\o,k 



f) 



(2nf^ 

Similarly, we find for the Hartree diagram 
\2g\Mm\ 



2u 2 + 



u 



— 3u arctan u 



u 



Ul?(p,kj 



dy 



s—xy 



d£ 



y{xy + " 



x[l + (xy + 



212 



2< + 



(u 2 - e) 



in 



(12) 



(13) 



ut\o,kf 



9\Mml 

(27TA) 4 



di 



2 , o<-3 



+(5u + 7£ + + 3£ 



1 + 



-u(7 + 3£ 2 ) 



In 



and for the Hartree diagram (d): 
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(15) 



with the auxiliary function a = t^y + \fu 2 — £ 2 + £ 2 y 2 . The limit x — > of eq.(15) is quite 
non-trivial since it requires a careful treatment of the singular function Ke(y + i0)~ 2 . The 
following representation is manifestly free of singularities, 

g\Mm\ /•« ^ f l dy f 3a + 2a 3 o _ 2£y(u 2 - £ 2 ^ 2 



(2tt/^ Jn 
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l y 2 I 1 + a 2 
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Note also the similarity of the expressions eqs.(ll,13,15) with the three-body single-particle 
potential given in eq.(ll) of ref.||. 




Fig. 3: The iterated In -exchange Fock graphs. Their isospin factor in symmetric nuclear matter 
is —3. 

Next, we come to the evaluation of the four iterated l7r-exchange Fock diagrams shown 
in Fig. 3. We start with diagram (e) with one medium insertion. Let us first consider the 
7riV-interaction vertices. The spin-dependent part of the matrix-product a ■ a a - b a ■ a a - b is 
Ha-b a- (ax b). In the case of graph (e) one makes the assignment: a = l + Q and b = I, with I 
the loop momentum and Q = Pi — p, where p belongs to the external nucleon line and pi to the 
internal nucleon line carrying the medium insertion. The important proportionality factor q is 
now not produced by the 7riV-interaction vertices as it was the case for the Hartree diagrams 
in Fig. 2. In order to isolate the factor q we combine the previously mentioned a ■ b term 
with the nucleon energy denominator (which in the actual calculation results from the dlo-loop 
integration) and employ the identity: I • (l + Q) [I ■ (l + Q + q = 1 — l-q [I - (l + Q + q The 
term coming along with 1 in this decomposition finally loop-integrates to zero. From the second 
term one can now easily isolate the factor q and take the limit to homogeneous nuclear matter 
of all remaining factors of the diagram. Putting all pieces together we end up with the following 
representation for the nuclear spin-orbit strength generated by the iterated l7r-exchange Fock 
diagram (e): 
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Note also that Uif(p,kf) in eq.(17) originates from the real part of the iterated l7r-exchange 
spin-orbit NN-amplitude (eq.(33) in ref. 0) evaluated in backward direction and integrated over 
a Fermi- sphere of radius kf. 

We continue with the computation of the Fock diagrams with two medium insertion, labeled 
(f) and (g) in Fig. 3. Diagram (f) with a symmetrical arrangement of the two medium insertions 



u 



2 + u 2 



u 2 

-(4 
2 v 



u 



u ) arctan 2(1 + u ) arctanw 



8 



leads to the following contribution to the nuclear spin-orbit strength: 
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with the auxiliary function t = xz + V^ 2 — x 2 + x 2 z 2 . At zero nucleon momentum (p = 0) the 
lengthy expression eq.(19) simplifies drastically to 
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Note the leading p 2 -behavior of the expressions in eqs. (12,20) derived from diagrams with two 
medium insertions. The results eqs. (8,10,18) belonging to diagrams with one medium insertion, 
on the other hand, show a leading linear dependence on the density p ~ u 3 . 

Finally, we have to evaluate the last two topologically distinct Fock diagrams in Fig. 3. 
Since they contribute equally to the nuclear spin-orbit strength we have given both diagrams 
the same label (g). In order to avoid very lengthy formulas we split their contribution to 
Ui s (p,kf) into a " factorizable" (g 1 ) and a " non-factorizable" part (g"). Technically these two 
pieces are distinguished by the feature whether the nucleon propagator in the denominator can 
be canceled by terms from the product of 7riV-interaction vertices in the numerator, or not. We 
find from the iterated l7r-exchange Fock diagrams (g) with two medium insertions the following 
"factorizable" contribution to the nuclear spin-orbit strength: 
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+ h x * + M V + 3x 2 _ 3u 2 _ 2) _ 3 (1 + ^ 2)2(1 + ^ 

2 24 



(21) 



which turns at zero nucleon momentum (p = 0) into the form: 

2g\Mm\ r £ 2 



zn VI TH C u r ( 

Ut\0,kf) = jtf^L *(i^^{«(3« a -9e 2 -17)+4(3 + e 2 )[arctaii(« + 



+ e) 2 

+ arctan(w - f )1 + (9£ 4 - 6w 2 £ 2 + 18£ 2 - 3m 4 - 26m 2 - 7)L(£, u) 



(22) 
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The analytical evaluation of the " non-factorizable" parts from diagrams (g) terminates with 
two non-elementary integrations and we find the following representation for their contribution 
to the nuclear spin-orbit strength: 



(47r/ ff )4 Jo * x 3 J-l 



a — arctan <r 



x - L(£,x) 



+ 
+ 



1 5 
- ln(l + a 2 ) + 4£y arctan a + 2u 2 - 2£ 2 - -a 2 



In 



\x + £y\ 



\x-£y\ 



4£y arctan cr + -(1 + £ 2 - x 2 ) ln(l + a 2 ) + 2u 2 - 2£ 2 



a 



1 \xR + {x 2 -i-e)iy\ \ 

R 11 \xR+ (1 + £ 2 -x 2 )£y\ J 



(23) 



with the auxiliary function i? = y(l + x 2 — £ 2 ) 2 + 4£ 2 (1 — y 2 ). At zero nucleon momentum 
(p = 0) eq.(23) turns into the (singularity free) form: 

Ag\Mml f u 



uL 9 "\o,k f ) 



(4ttA) 4 Jo 



'! + ?)* U y 2 



dy 



^l2Z(l + e-2y 2 



arctan a — a 



ln(l + a 2 ) - ln(l + u 2 - £ 2 ) - 2£ya + 



+y 



+2f (1 + C )(l + y ) Ju 2 - £ 2 - arctan Ju 2 - £ 2 



2e(i+e)(u 2 -z 



l + u 2 -£ 2 



(24) 



Note again the similarity of the expressions in eqs. (19,21,23) with the three-body potential 
given in eq.(13) of ref.||. We also like to emphasize that the techniques used to reduce six- 
dimensional principal value integrals over the product of two Fermi-spheres of radius kf to 
at most double integrals have been checked rigorously in ref.[§] via the Hugenholtz- van-Hove 
theorem. 

Let us end this section with power counting considerations for the spin-dependent inter- 
action energy S spin . Counting the quantities p, q, kf and m n collectively as small momenta, 
we deduce that the relativistic l/M 2 -correction from the l7r-exchange Fock graph is of fifth 
order, while all contributions from iterated l7r-exchange are of fourth order in small momenta. 
Irreducible 27r-exchange gives also rise to a fifth order contribution to the (real part of the) 
single-particle potential (see eq.(16) in ref.|§). To the spin-orbit NN- amplitude irreducible 2n- 
exchange contributes as a (higher-order) relativistic 1/M-correction (see eqs. (22,23) in ref.||) 
and it will therefore enter E sp j n at sixth order in small momenta. We have also checked that the 
relativistic 1/M-corrections to iterated l7r-exchange start to contribute to S spin first at sixth 
order in small momenta. From all that we can conclude that the present calculation of the 
nuclear spin-orbit strength Ui s (p, kf) is complete up-to-and-including third order in small mo- 
menta. Strictly speaking, this counting argument applies only to the long-range effects induced 
by chiral In- and 27r-exchange since it does not cover the possible short-range contribution, 
Uig hort \p, kf) = Cisp, which is of course also of third order in small momenta. From this point 
of view the oxu-exchange term eq.(5) (together with the Fock contributions) provides a model 
for the a priori undetermined low-energy constant C/ s . In order to learn about the convergence 
of the chiral expansion of the nuclear spin-orbit strength Ui s (p,kf) one should calculate the 
contributions from irreducible 27r-exchange of fourth order in small momenta. Work along this 
line is in progress. Furthermore, we note that spin-orbit strength generated by the three-body 



force diagrams in Fig. 2 of ref . |l7j is (formally) of higher order in small momenta. These should 
also be evaluated with the formalism introduced in section 2. 
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Diagram 


l7r-Fock 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g') (g") 


U ls (0,k fo ) 


0.40 


-52.96 


37.35 


41.11 


12.57 


30.33 


14.69 


-44.50 -3.85 


Uis(k f0 , kfo) 


-0.28 


-40.05 


31.63 


13.95 


-35.35 


25.63 


4.35 


-29.05 8.04 



Table 1: Contributions of individual diagrams to the nuclear spin-orbit strength Ui s (p,kf ) 
at p = and at p = km = 272.7 Me V. The units are MeVfm 2 . 



3.1 Results 

For the numerical evaluation of the nuclear spin-orbit strength Ui s (p, kf) we use consistently the 
same parameters as in our previous works |7|, §. We choose the value g& = 1-3 for the nucleon 
axial vector coupling constant. The weak pion decay constant has the value f n = 92.4 MeV 
and M = 939 MeV and m n = 135 MeV are the masses of the nucleon and the (neutral) pion, 
respectively. 

In the second row of Table 1, we present numerical values for the contributions of indi- 
vidual diagrams to the spin-orbit strength Ui s (0,kfo) at nuclear matter saturation density 
kfo = 272.7 MeV. As expected the relativistic 1 / M 2 -correction from the l7r-exchange Fock graph 
is a very small 1.1% effect. The contributions of individual iterated l7r-exchange diagrams are 
surprisingly large. In several cases they even exceed the empirical value Ui s ~ 35 MeVfm 2 in 
magnitude and moreover they are of alternating signs. The basic reason for these large values 
is the large scale enhancement factor M (the nucleon mass) entering the iterated l7r-exchange. 
The proportionality factor M stems from the energy denominator of such second-order dia- 
grams which is a difference of small nucleon kinetic energies. Adding up the entries in the 
second row of Table 1 one gets Ui a (0,kf ) = 35.1 MeVfm 2 , which is in perfect agreement with 
the empirical value of the nuclear spin-orbit strength Ui s ~ 35 MeVfm 2 JTl|, [12]]. The predicted 
total sum is dominated by the contribution L^ s (0, kfo) = 38.1 MeVfm 2 of the four Hartree 
diagrams (a), (b), (c) and (d) (see Fig. 2). Interestingly, the same feature, namely the numeri- 
cal suppression of the iterated l7r-exchange Fock diagrams against the Hartree diagrams, holds 
also for the (spin-independent) nuclear mean-field U{p,kf) studied in ref. |5J. The novel spin- 
orbit strength generated (almost completely) by iterated l7r-exchange is neither of relativistic 
nor of short range origin. It is in fact linearly proportional to the nucleon mass M and its 
inherent range is the pion Compton wavelength m~ x = 1.46 fm. The latter feature tempts to 
an unconventional interpretation of the strong nuclear spin-orbit interaction. The potential 
Vis underlying the empirical spin-orbit strength C7j s = VJ S rf s ~ 35 MeVfm 2 becomes a rather 
weak one, namely Vi s ~ 17 MeV, after the identification of the effective range r/ s with the pion 
Compton wavelength, r/ s = m" 1 = 1.46 fm, as suggested by the present calculation. 

As a side remark we consider Ui s (0, kf) in the chiral limit m n = 0. In that case all occurring 
integrals can be performed analytically and we find the following simple expression: 

g\Mk 2 f /4 tt 2 \ g 2 A k 3 f 
U ls (0, k f ) ]nhr=0 = - - - + QfA l A f f — 2 , (25) 



(4tt^) 4 V3 2 J 3(4ttAM) 2 

which gives at kf = 272.7 MeV the negative value —15.0 MeVfm 2 . 

In Fig. 4, we show the pion mass dependence of the nuclear spin-orbit strength Ui s (0, fc/o) at 
nuclear matter saturation density (and zero nucleon momentum). One observes a sign change 
and a pronounced maximum at m n ~ 100 MeV. 
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Fig. 4: The pion mass dependence of the nuclear spin-orbit strength Ui s (0,kfo) at saturation 
density k f0 = 272.7 Me V. 
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Fig. 5: The pion mass dependence of the depth of the spin-independent single-particle potential 
U(0,kf ) j^j at saturation density kf = 272.7 MeV. 
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This strong m^-dependence of Ui s (0, kfo) has its origin in the alternating signs of the large 
contributions from individual diagrams as well as their own specific m^-dependence. For com- 
parison, we show in Fig. 5 the pion mass dependence of the depth of the (real) single-particle 
potential U(0,kf ) calculated in ref. ||. One observes a weak variation of the potential depth 
U (0, kfo) by at most 10% if m n runs from zero (chiral limit) to the physical value m n = 135 MeV. 
With increasing pion mass the attractive nuclear mean-field generated by chiral 27r-exchange 
gets however soon lost and it turns into repulsion above ~ 215 MeV. 

In Fig. 6, we show by the full line the dependence of the calculated nuclear spin-orbit strength 
Ui s (0,kf) on the nucleon density p = 2k^/3ir 2 . One observes in the region p < 0.4 fm~ 3 an 
approximate linear growth of Ui s (0,kf) with the density as it is, for example, known from o~- 
and a;-exchange (see eq.(5)). 

In a finite nucleus the spin-orbit force acts mainly on the surface where the density gradients 
are largest and the density has dropped to about half of the central density. The replacement 
f( r ) V/(r) — > | v /(r) (instead of f(r) V/(r) — > V/(r) valid for weakly inhomogeneous nuclear 
matter) describes then more realistically the situation for a finite nucleus. The dashed line in 
Fig. 6 shows the spin-orbit strength which results if the contributions from the diagrams with 
two medium insertion (b), (c), (d), (f) and (g) are weighted with a factor 1/2. This different 
weighting leads to a substantial reduction of the total spin-orbit strength such that only about 
18% of the empirical value U\ B ~ 35 MeVfm 2 are left at nuclear matter saturation density 
p = 0.178 fm- 3 . 




-10 

0.0 0.1 0.2 0.3 0.4 

P [fm 3 ] 

Fig. 6: The full line shows the nuclear spin-orbit strength Ui s (0, kf) at zero nucleon momentum 
(p = 0) versus the nucleon density p = 2/c 3 /37r 2 . If the diagrams with two medium insertions 
are weighted with a factor 1/2 the dashed line results. 

In Fig. 7, we show by the full line the dependence of the calculated nuclear spin-orbit 
strength Ui s (p,kfo) at saturation density kfo = 272.7 MeV on the nucleon momentum p for 
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< p < kjQ. One observes a very strong p-dependence which leads even a sign change of 
U[ s (p, kfo) above p = 200 MeV. Again, this strong p-dependence of Ui s (p, kf ) has its origin in 
the alternating signs of the large contributions from individual diagrams as well as their own 
specific p-dependence. The numerical values in the second and third row of Table 1 indicate how 
these contributions from individual diagrams change with the nucleon momentum from p = 
to p — kfo- One also should keep in mind that the scale relevant for momentum dependences 
is here set by the pion mass, m n = 135 MeV. Note that the ratio p/m n changes by two units 
from p = to p = kf . Furthermore, the dashed line in Fig. 7 corresponds to the weighting of 
diagrams with two medium insertions with a factor 1/2. It shows the same strong p-dependence 
including a sign change as the full line in Fig. 7. 



40 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




p [MeV] 

Fig. 7: The momentum dependence of the nuclear spin-orbit strength Ui s (p,kfo) at saturation 
density kfo = 272.7 MeV. The full (dashed) line corresponds to the weighting of diagrams with 
two medium insertions with a factor 1 (1/2). 

Let us briefly summarize our results. We have demonstrated here that the nuclear spin-orbit 
interaction is not necessarily a relativistic effect. Iterated l7r-exchange (i.e. loosely speaking, the 
l7r-exchange spin-spin and tensor force in second order) generates large nuclear spin-orbit terms 
which in fact scale linearly with the nucleon mass M. For weakly inhomogeneous nuclear matter 
and small nucleon momenta (p < 70 MeV) the spin-orbit strength from iterated l7r-exchange 
agrees (at saturation density) well with the empirical shell model value Ui s ~ 35MeVfm 2 . The 
strong p-dependence of Ui s (p, fc/o) (including a sign change) and the different weighting of V/(r) 
and f(r) V/(r) at the nuclear surface leave however questions about the ultimate relevance of 
this 27r-exchange spin-orbit interaction for real nuclear structure. Nuclear structure calculation 
which use the calculated spin-orbit strength Ui s (p, kf) as input are necessary in order to clarify 
the role of the spin-orbit interaction generated by 27r-exchange. 
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4 Isovector single-particle potential 



In this section we generalize our previous calculation || of the single-particle potential to isospin 
asymmetric (homogeneous) nuclear matter. Any relative excess of neutrons over protons in 
the nuclear medium leads to a different single-particle potential for a proton and a neutron. 
This fact is expressed by the following decomposition of the single-particle potential in isospin 
asymmetric nuclear matter: 



U(p,k f )+iW{ P ,k f )- UrfakA+iWfakf) r 3 S + 0(S 



(26) 



Here, U(p, kf)+i W(p, kf) is the (complex) single-particle potential in isospin symmetric nuclear 
matter. The term linear in the isospin asymmetry parameter 5 = (N — Z)/ (N + Z) defines the 
(complex) isovector single-particle potential Ui(p, kf) + iWi(p, kf), and r 3 — > ±1 for a proton 
or a neutron. The situation of an isospin asymmetric nuclear medium is realized by the simple 
substitution: 

d(kf-\Pj\) - l -^e{k p -\p j \) + l -^9{k n -\p j \), (27) 

in the medium insertion eq.(4). Here, k Pjn = kf(l =)=5) 1 / 3 denote the (different) Fermi momenta 
of protons and neutrons. Differences in comparison to the calculation of U(p, kf) + i W(p, kf) 
in ref . || occur only with respect to isospin factors and the radii of the Fermi-spheres, k p ^ n = 
^/•(l^^) 1 / 3 . In practise the isovector single-particle potential Ui(p, kf) +i Wi{p, kf) is obtained 
by differentiating the r 3 -components of the diagrammatic expressions with respect to 5 at 5 = 0. 



4.1 Real part 

Without going into further technical details we enumerate now the individual contributions to 
the real part of the isovector single-particle potential Ui(p, kf). 

i) l7r-exchange Fock diagram in Fig. 1 including the relativistic l/M 2 -correction: 

,2 



U^\p,k 



f) 



2L(x,u) 



1 3 2 

3(4ttA) 2 

-(u 2 + x 2 )L(x,u) 4 



77T 

2M + # 



U[U 



u[u — X 



X 
2\2 



2[1 + (u + x) 2 }[l + (u-x) 2 ] 
with x = p/m n and u = kf/m n . The auxiliary function L(x,u) has been defined in eq.(7). 
ii) Iterated l7r-exchange Hartree graphs in Fig. 2: 

g\Mm 4 K u 2 ( /u _\ /u ^\ , 5 



(28) 



ui a) (p,k f ) 



\( l) arctan(u — x) — f — + l) arctan(w + x) H — L(x,u)\ 

487^/^ l v x ' v x ' 2 J 

3(47T/ 7r ) 4 J-i I 



(29) 



U?\p,k 



f) 



2s 2 + s 4 



s 5 s' 



1 + 8 



2\2 



u — xy 
2uxy + {u 2 — x 2 y 2 ) In 



1 + s 2 
u + xy 



21n(l + s 2 



u — xy 



(30) 



with s' = uds/du and the auxiliary function s has been defined after eq.(ll) 



U<t c \p,kj 



3(47^ 
2s' 



dy 



s-xy 



(xy + E 



s [i + (^ + 2 ] 2 



s'(2£ + xy) 



ein 



- 2u 



(xy + 2 



2u£ + {u 2 - £ 2 ) In 



u+J, 



- (2s + s')u 2 In 



2 ,._ w + £ 



(31) 



15 



Ul d \p,k f ) = 



3(4ttA) 4 

J2 



Jo X J-l I 



2x£y 



x 2 — £ 2 y 2 



TT ^-21n(l + ^) 



+ 



■(5cr 4 + 9a 2 + 6 - 4crV) - 6 ln(l + a 



(3 2 ) 



.(l + a 2 ) 2 

with a' = uda/du and the auxiliary function a has been defined after eq.(15). The symbol 
■f-i dy stands for a principal value integral. 

iii) Iterated l7r-exchange Fock graphs in Fig. 3: 

5g\Mm^u 2 



ui e \ P ,k f ) = 



6(4vr)3/4 



2u 



(u+x)/2 



(u-x)/2 X(l + 2£ 2 ) 



X 



(1 + 4£ 2 ) arctan2£ - (1 + 8£ 2 + 8£ 4 ) arctan£ 



(33) 



u { /\ P ,k f ) 



g\Mm\ ( uG(x,u) dG(x,u) 



3(47^) 



4x 2 



du 



+ / dy I dz 



1 yz9(y 2 + z 2 - 1) sV 



t 2 -ln(l + t 2 



-l J-i \yz\\Jy 2 + z 2 — 1 1 + s 2 
with the auxiliary function G(a;, u) = u(l + u 2 + x 2 ) — [1 + (u + x) 2 ][l + (u — x) 2 ] L(x, u) 



(34) 



U ( /\p,k f ) = ff^l { G(u, u) (1 + x 2 - u 2 )L(x, u)-u 



3(47r/ w 

ru 

+5u / di 
Jo 



:(l + X 2 -e)L(x,0-l 



du 



}■ 



(35) 



uP(p,k,) 



g\Mm\ /■« £ 2 r 1 



(4tt^) 4 7o 



6X J-l 



a 



l + a 



2 (3 + 5a 2 + 8(7(7') - 3 ln(l + a 2 ) 



x 



X 



ln k + + 1 ln 1^ + (^ 2 - e 2 - i)^l 



k -ez/i \xR+(e + i-x 2 )ty\ 



£ ln(l + cr 2 ) - a' 



R! \xR+ (x 2 -£ 2 - , - x 2 + 3£ 2 ) - xR' 



R 2 \xR + (£ 2 + 1 -x 2 )£y\ R[xR + (x 2 - £ 2 - l)£y] 
y(l -x 2 + 3£ 2 ) +xR' 2xy 



R[xR + (£ 2 + 1 - x 2 )iy] x 2 -i 2 y 2 . 

with R! = dR/dl; and the auxiliary function R has been defined after eq.(23). 
iv) Irreducible 27r-exchange Hartree and Fock graphs: 



(36) 



(37) 



j(0 = 3(1+2^+ 5^i) l^^+yi+e) 



+2 



5 + 26s 2 - 79g\ + 2e 2 (i + Wg 2 A - Wg A A )]ZyJT+?MZ + ^ + 



+ (17 + 242«? 2 - 787g\)e ~ (3 + 14<& + 15sl)£ 4 



+ 



60(1 + 6g 2 A - 15g\)e + 4(1 + Wg 2 A - 59g A A )e 



ln 



2A 



(38) 
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v) Power divergences specific for cut-off regularization: 



ui A \p,kj 



2A fcj 



26g\M + 5(3g 2 A + l)(l-g A )A 



(39) 



The term linear linear in the cut-off A stems from iterated l7r-exchange with a contribution 
of the Hartree diagram (a) and the Fock diagram (e) in the ratio 8:5. The term quadratic 
in the cut-off A, on the other hand, originates from irreducible 27r-exchange. Note that the p- 
independent contribution to Ui(p, kf) in eq.(39) is just twice its contribution to the asymmetry 
energy A(kf) (see eq.(29) in ref. 0]). This relative factor of 2 is typical for a momentum 
independent NN-contact interaction, to which the power divergences are completely equivalent, 
as emphasized in ref.pj. We use consistently the value A = 646.3 MeV ~ 7f^ of the cut-off scale 
which has been fine-tuned in ref.[[F) to the binding energy per particle, —E(kfo) = 1 5. 26 MeV. 




Fig. 8: The real part of the isovector single-particle potential Ui(p,kf) in isospin asymmetric 
nuclear matter versus the nucleon density p = 2k^/3TT 2 . The full (dashed) line corresponds to 
p = (p = kj). 

In Fig. 8, we show by the full line the total (real) isovector single particle potential {7/(0, kf) 
of our calculation for a nucleon at rest (p = 0) as a function of the nucleon density p = 2kj/37r 2 . 
The shape of this curve is very similar to the asymmetry energy A{kf) (see Fig. 7 in ref.0). In 
comparison to A{kf) the scale on the ordinate is stretched by a factor of about 1.4. Interestingly, 
the (real part of the) isovector single-particle potential {7/(0, kf) has its maximum close to the 
saturation density po = 0.178 fm~ 3 . The actual value at that point is {7/(0, kf) = 47.0 MeV. 



This prediction is comparable to the value U\ ~ 33 MeV |TT[] used in shell model calculations 
or the value Ui ~ 40 MeV |13| deduced from nucleon- nucleus scattering in the framework 
of the optical model. The dashed line in Fig. 8 shows the density dependence of the (real) 
isovector single-particle potential Uj(kf,kf) at the Fermi surface p = kf. At that point the 
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(real) isovector single-particle potential comes out always more repulsive than at p — 0. Note 
also that the (possibly unrealistic) downward bending branches of the curves in Fig. 8 start at 
densities higher than those relevant for conventional nuclear physics. 

Furthermore, we show in Fig. 9 the momentum dependence of the (real) isovector single- 
particle potential Uj(p,kf ) at saturation density kfo = 272.7 MeV. The p-dependence of 
Ui(p,kfo) is non-monotonic in the interval < p < kfo- One observes a broad maximum 
at p = 230 MeV where Ui(p, kfo) has increased by about 30% to the value 62.9 MeV. Note that 
in comparison to the spin-orbit strength Ui s (p, kf ) shown in Fig. 7 the p-dependence of the real 
part of the isovector single-particle potential Ui(p, kf ) is very moderate. Most useful would be 
nuclear structure calculations using the calculated isovector single particle potential Ui(p, kf) 
as input. 
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Fig. 9: The momentum dependence of the real part of the isovector single-particle potential 
Ui(p,kf ) at saturation density kf = 272.7 MeV. 



4.2 Imaginary part 

In this subsection, we discuss the imaginary part Wi(p, kf) of the isovector single-particle 
potential. According to eq.(26), it splits and shifts the half-width of a neutron-hole state and a 
proton-hole state in the Fermi-sea (with momentum < p < kf) by the amount ±5 Wi(p, kf). 
Within the present calculation the imaginary part Wi(p,kf) arises entirely from iterated one- 
pion exchange. It is advantageous to present analytical formulas only for complete classes of 
diagrams. We find from the iterated l7r-exchange Hartree diagrams in Fig. 2 (including the 
analogous graph with three medium insertions) the following contribution to the imaginary 
isovector single-particle potential: 

wi H \p,k f ) = f|^|2 M 4 -10 M 2 -y M V + 8 M 2 ln(l + 4 a ; 2 ) 
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(40) 



The iterated l7r-exchange Fock diagrams in Fig. 3 (again including the analogous graph with 
three medium insertions) lead, on the other hand, to the following expression: 
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(41) 



with the auxiliary function = £y + -\/w 2 — re 2 + £ 2 y 2 . If the delta-function 6(x — £\y\) is 
used to eliminate the rfy-integration in eqs. (40,41) the remaining d£-integral extends over the 
restricted region x < £ < u. The sum of both contibutions eqs. (40,41) evaluated at zero nucleon 
momentum (p = 0) can even be written as a closed form expression: 



W T (0, kj 



g\Mm\u 2 f 9m 6 + 40m 4 + 27m 2 12m 4 + 40m 2 + 27 



384tt 3 /4 



2(1 + M 2 ) 2 
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(42) 



Finally, we show in Fig. 10 the momentum dependence of the imaginary isovector single- 
particle potential Wi(p, kfo) at saturation density kfo = 272.7 MeV. The associated value at zero 
nucleon momentum, W/(0, kfo) = 26.7 MeV, agrees within 10% with the isoscalar half-width 
W(0, kfo) = 29.7 MeV found in ref.[§. As a consequence of the decreasing phase space available 
for redistribution of a nucleon-hole state's energy, the curve in Fig. 10 drops with momentum p 
and Wi(p, kfo) reaches zero at the Fermi-surface p = kfo. The exact vanishing of Wi(p, kf) at 
the Fermi-surface p = kf is even separately true for the class of iterated l7r-exchange Hartree 
diagrams and the class of iterated l7r-exchange Fock diagrams. The conditions Wj (kf, kf) = 
and Wf'(kf,k f ) = serve as an excellent (analytical and numerical) check on the involved 
calculations leading to eqs. (40,41). 

In summary, we find that the predictions from chiral In- and 27r-exchange for the real part 
of the isovector single-particle potential Uj(p,kf) agree fairly well with empirical values. The 



calculated imaginary part Wi(p, kf) fulfills the constraints imposed by Luttinger's theorem [|14 
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Fig. 10: The momentum dependence of the imaginary part of the isovector single-particle po- 
tential Wi(p, kfo) at saturation density kfQ = 272.7 MeV. 
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